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PROPULSION ESTIMATION TECHNIQUES 


by 

Enoch C. Temple 

Associate Professor of Mathematics 
Alabama A&M University 
Normal, Alabama 

ABSTRACT 

This paper gives a brief overview of statistical tools that are 
needed to perform post flight /test reconstruction of state variables. 
Linear regression, recursive linear regression, and the exact connection 
between the Kalman filter and linear regression are discussed. The 
regression connection is expected to serve as an aid in the application 
of a recently developed analytical method of flight reconstruction to 
single engine test firing data. 
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1. INTRODUCTION 


This paper considers the issue of post flight /test state variable 
reconstruction through the application of observations made on the 
output of the Space Shuttle propulsion system. The incentive for the 
concepts discussed here is rooted in a NASA contract report by Rogers 
(1987). Rogers used the Kalman filtering procedure to reconstruct the 
state variables of the Space Shuttle propulsion system. 

Since its development by Kalman (1960), the Kalman filtering metho- 
dology has been a popular procedure with control engineers. Within recent 
years , it has seen increased applications in other engineering areas as well 
as in nonengineering areas. Although the applications have been diverse 
and research on some of its defects has been quite extensive, applied statis- 
ticians have not taken full advantage of the tool. This lack of use by some 
applications areas may be because many publications on the topic have 
obscured statistical simplicity . Often times the basic statistics concepts are 
embedded in a discussion of numerical solutions to differential equations . 

An objective of this paper is to give the general setup of the 
Kalman filter and its connection to linear regression. A second, objective 
is to examine the Rogers (1987) reconstruction methodology for applica- 
tion to the reconstruction of the state vector of a single Space Shuttle 
Main Engine (SSME) by using static test firing data. 

Throughout this paper, underlined capital letters are used to 
denote vectors, capital letters denote matrices, and the identity matrix 
is denoted by the letter I. All vectors are of the column type and the 
transpose of any matrix or vector is denoted by using the letter T at the 
superscript position. The letter E denotes the expectation operator, 

N(U, Z) denotes the multivariate normal probability distribution with 
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mean vector U and covariance matrix Z. The caret symbol written 
directly above a scalar or vector denotes a statistical estimator of that 
vector or scalar. 

2. ESTABLISHMENT OF NOTATION 

For any Space Shuttle flight, let denote the observed values of 
vector Z at time t. Each component of Z^ represents a relevant output 
of the Space Shuttle propulsion system that can be measured. For 
example, the components of Z^. may be chamber pressure, oxygen flow 
rate, hydrogen flow rate, etc., for each of the three SSME’s. Rogers 
(1987) lists 71 components for vector Z^ and 35 components for vector 
Xj . where Xj. is defined below. 

Vector X^ is a state vector of parameters to be estimated at time t. 
It is assumed that the observation vector Z^ is a function of the state 
vector X.. That is. 


5t = 


cx t ) 


+ v. 


( 2 . 1 ) 


where h is some function and V t N(0^,R t ). State vector X^ is known 
to change with respect to time according to the equation 


x t = f (X, 


,t) + w. 


( 2 . 2 ) 


where f is some function, W t v N(0^,Q t ) and ) = 0. 

If it is assumed that f and h in equations (2.1) and (2.2) are 
linear, then numerical procedures allow us to transform these equations 
to the form 


— K “ F K £k-1 + 


(2.3) 
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” H K — K + — K 


(2.4) 


where K represents discrete values of t, W*. v N(0, Qt.) and V„ v 

—iV — iv — 

N(0,R k ). 

After a Space Shuttle flight has taken place, the post flight reconstruc- 
tion procedure seeks to use the observed values of and the equations in 
(2.3) and (2.4) to reconstruct (estimate) the parameters X R so that the estima- 
tion error is minimized. The estimated value of X T , is denoted bv X T r. 

3. A RECURSIVE ESTIMATOR OR X 

In order to complete the development of the reconstruction process 
mentioned in section 2, we consider a system of equations 

Z = H X + e (3.1) 

where (3.1) is a system of linear equations that has been generated by 
making n observations on the single equation 

Z* = H* X + e* (3.2) 

where Z* and e* are vectors or scalars. H* is a jxp matrix where 
j = number of components in vector Z* and X is fixed. Vector X 
being fixed means that X does not change with time or does not change 
as the number of observations on Z increase. If Cov( e) = Z and 
det ( Z ) 4 0 , then the weighted least squares estimate of vector X is 
given by 

X = (H T z" 1 H) -1 H T z" 1 Z (3.3) 

and the covariance of the estimate is given by 
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Cov(X) = E [(X - X) (X - X) T ] = (H T E _1 H)" 1 = P . (3.4) 

The reader should be reminded that with the appropriate assump- 
tions, the maximum likelihood, the minimum variance unbiased linear as 
well as the minimum mean- squared error estimator of X are identical to 

A 

the estimator given by equation (3.3). We also remark that Cov(X) is a 

A 

measure of the quality of the estimates in vector X. A detailed dis- 
cussion of these estimators may be found in Elbert (1984). 

In the sequel that follows, the estimator in equation (3.3) will be 
rearranged so that X can be estimated by using a recursive process. 

The development of the recursive process is started by assuming that 
there are K observations on equation (3.2) which form the system 

Zj =H K X + £ K (3.5) 

where Cov = X K and det (Ir) 4- 0. By equation (3.3) 

Xr = ( h k T ^k 1 h k } 1 h k z k X hi ' (3 ‘ 6) 

is the estimate of X where observations up to and including observation 
K are used. The covariance of the estimator is denoted by and 

P K = E[(X r - X) (>k - X) T ] = < Ir 1 H k ) _1 . (3.7) 

Assume that an additional measurement hi+1 has been made. Let 

^K+l = H K+1 - + ^K+l (3,8) 

T T 

where Cov (jEr +1 ) = z k+ 1 311(1 E -K+l^ = 0> Note that E (iK -^C+l^ 
means that the (K+l)th observation is independent of the 1st K observa- 
tions. Putting the K observations with the (K+l)th observation yields 
the system 
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I I - 


\z 


— K+l 



H 


K+l 


X + - 


Vs 


rK+1 


(3.9) 


where 


/% \ 

Cov ( I 

Vk+i/ 


J K 


J K+1 


) 


By substituting into equation (3.3) we get 


^K+l 


(H 


K 


, -i ; 

'K i 


H K+1 ) » 


0 V h k Yl" 1 

iAh Zi 


J K+1 K+l 


Note that we can write 


-K+l " + ( -K+l " ^K ) * 


(H k ! h Li> 


/ Z K < 0 V^K \ 

' 0 ^k+Ak+i' 


(3.10) 


(3.11) 


Substituting results from (3.6) and (3.10) into (3.11) and using matrix 
algebra, we get 


^K+l ^K + (P K X + H K+1 H K+1 ) 1 H K+1 Z K+1 ( -K+l " H K+1 ^K ) * (3 * 12) 

By substituting into equation (3.7) and using matrix algebra, we get 


K+l 


T > T 

< h k : W 


( z 'k i o \/ h k VT 
' 0 ! e k+ih 


(P K + H K+1 £ K+1 nR+l* 


K+l 


Substituting P R+1 into equation (3.12) yields 
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(3.13) 


“ ?K + P K+1 H K+1 Z K+1 ( -K+l H K+1 ^K* 

and 

P K+1 = ^R 1 + H K+1 Z K+1 H K+1 ) 

Equations (3.13) and (3.14) provide recursive equations for estimating 
X and its corresponding covariance matrix P. 

Recall that the estimator given in equation (3.3) may be identified 

/K 

by using the symbol X^. Therefore, if X is fixed as it is in equation 
(3.1), the recursive equations of (3.13) and (3.14) have no statistical 
advantage over the estimator given in equation (3.3). However, if X 
changes with time where there is one observation per time interval, 
then the recursive equations are quite useful. In fact, equations (3.13) 
and (3.14) when combined with equations (2.3) and (2.4) form the 
Kalman filtering process. 


4. THE KALMAN FILTER 


It is worthwhile to mention that Brown (1983) and Gelb (1974) 
used matrix theory to express equations (3.13) and (3.14) as 


P K+1 " (I " K K+1 H K+1 ) P K 

^K+l ^K + K K+1 ( -K+l ' H K+1 ^K J 
where 


K K+1 “ P K H K+1 (Z K+1 + H K+1 P K H K+1 ) 


and K k+1 is called the gain matrix. 


(4.1) 

(4.2) 


(4.3) 
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For the convenience of the reader, equations (2.3) and (2.4) are 


reprinted here. That is. 


— K+l F K — K + 


W 


K 


(4.4) 


and 


-K+l H K+1 -K+l + -K+l * 


(4.5) 


When equations (4.1) and (4.2) are combined with equations (4.4) and 
(4.5), two types of estimates of X R+1 are possible for each K. A 
notation by Gelb (1974) allows for the two estimators to be distinguished. 
That notation is 

X R+1 7 = the estimate of X R+ ^ using all observations up to and 
including observation K. 

P R ^.^ = the estimate of P using all observations up to and including 
observation K. 

V 

a (+) 

X K ^' = the estimate of X R+ ^ using all observations up to and 
including observation (K+l). 

The values of X R ^^ and P^^ are obtained by using equations (4.4) 
and (4.5). The X R+1 vector is often called the extrapolated or pre- 
dicted value of X R+1 and P K +j^ is called the extrapolated variance of 



A ( 'i A ( “ ) 

The values of X^^ and P R ^ } are computed by substituting Xr+I 
for X R in equation (4.2) and Pj^j* for P R in equation (4.1). 

The computation summary is 

?K« > = F K 2k +) < 4 ' 6 > 
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(4.7) 


( ^ = F„ P^, +) + 


K+l 


K K K 


Qi 


y ("*") — v ( ) + y (7 - H X ^ 

-K+l _ -K+l n K+l ( -K+l n K+l -K+l ; 


P K« = (I - K K+1 H K + 1> P Kn 


K K+1 ~ P K+1 > h k+i <h k+i p k+i > h k+i + R K+l’ 


(4.8) 

(4.9) 

(4.10) 


Equations (4.6) to (4.10) completely describe the Kalman filtering process 
when the original functions f and h of equations (2.1) and (2.2) are 
linear. 

However, when functions f^ and h are nonlinear, the state vectors 
X R (K = t Q , t v ... t T ) are estimated through the application of an 
extended Kalman filtering procedure. This procedure requires a linear- 

j|e 

ization of functions f and h about some known state value X R . The 
next paragraph provides a brief overview of the extended Kalman filter 
concept . 

* 

Let X„ be some known value of X and assume that AX is small. 

A first degree Taylor series approximation of functions (2.1) and (2.2) 
may be 


* * AX + W t (4.11) 

K 

and 

* • AX + V . (4.12) 

X=X K 

* . * * 

If it is assumed that X R is selected so that X R = f(X R , K), then 
equations (4.11) and (4.12) become 


Z t * h(X R ) + 



X R + AX * f(X*, K) + 


lax 


x 
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(4.13) 


AX 



AX + W K 


and 



(4.14) 


Equation (4.13) is called the linearized dynamics equation and (4.14) 
is the linearized measurement equation. Note that equations (4.13) and 

(4.14) may be transformed to equations that are equivalent to equations 
(2.3) and (2.4). Hence for each discrete time K, AX K can be estimated 
and error covariance matrices can be determined. The state vector 
estimate at time K is then given by 


X ( +) - x * + AX (+) 
-K -K -K 


A ( +) 

where AX^ ' is computed by substituting from equations (4.13) and 

a 

(4.14) into equation (4.8). Vector X^ + ^ is computed by letting 
* 

-K+l ~ -K and repeating the above procedure. 


The extended Kalman filter has performed well in a large class of 
applications. However, there are occasions when divergence occurs in 
the state vector estimates. Divergence occurs when the computed 
entries of the error covariance matrix P R become small as compared to 

the actual error in the estimate of the state vector. The cause of this 

/ 

divergence is not due to a defect in the filtering procedure, but may 
be caused by the linear approximation procedure, numerical rounding 
error, or an inadequate model of the system being studied. Additional 
possible causes of divergence are mentioned by Gelb ( 1974) . 
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A publication by Varhaegen and Van Dooren (1986) is representa- 
tive of the theoretical and experimental analyses that are currently 
being done on the divergence problem. The Rogers (1987) computer 
code has implemented the U-D factorized algorithm as a means of con- 
trolling numerical roundoff error that may lead to divergence in the 
state vector. Checks for other sources that may generate divergence 
are conducted through computer evaluations. 

5. CONCLUSION 

At this point we have reviewed the general setup for a regression 
problem where the parameters to be estimated are fixed. This lead to 
recursive equations (4.1), (4.2), and (4.3) which allowed for the 

x. 

development of an estimation procedure for a time varying parameter. 
When functions f and h are linear, it has been clearly stated that 
equations (2.3) and (2.4) are the essential ingredients for state 
variable reconstruction. If f and h are nonlinear, X R can be estimated 
by linearing h and f about some known vector and treating the 
linearized equations as if they are equations (2.3) and (2.4). 

The problem of applying the extended Kalman filtering procedure 
to the static test firing data remains. The basic approach for the 
application is identical to the presentation given in Section 4. There- 
fore, the computer codes that have been prepared by Rogers (1987), 
which are currently operational on the MSFC computer system, may be 
modified so that static test data may be analyzed. The static test 
setup differs from the actual flight data in that many flight associated 
modules will become inactive. Therefore, after adjustments are made 
for the analysis of static test data, it will be necessary to evaluate the 
performance of the computer code for the divergence of parameter 


estimates. 
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